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Abstract 

We show that the recently derived (g-) discrete form of the Painleve VI equation can be 
related to the discrete Pud i n particular if one uses the full freedom in the implementa- 
tion of the singularity confinement criterion. This observation is used here in order to 
derive the bilinear forms and the Schlesinger transformations of both q-Pm and g-Pyi- 



1. Introduction. 



In a recent paper [1], the first two authors obtained a first order system of (/-difference 
equations the continuous limit of which was the sixth Painleve equation. This completed 
the list of discrete forms of the Painleve transcendents, since a discrete Pyi was the only 
one that was missing. The method used for obtaining this mapping was to consider a q- 
difference version of the monodromy preserving deformation theory of linear differential 
equations. (We recall the strong connection of the continuous Painleve equations [2] to 
the latter approach.) The (/-difference system that was called g-Pyi is the following: 

a 3 a 4 (x - q n b 1 )(x - q n b 2 ) . s 

VV -= U /,.,)U- 6.) (Lla) 

xx = Wy-g"ai)(i/-g"a2) 
{y - a 3 ){y - a 4 ) 

where x = x(n), x = x(n + 1), x = x(n — 1) and a\, . . . , 64 are constants constrained by 
aia 2 /a 3 a 4 = (76162/6364. This system is integrable in the sense that it results from the 
compatibility of linear (/-difference equations. Moreover, the singularity confinement [3] 
property of (1.1) was verified explicitly. 

In [4], the remaining authors, while missing the connection with Pyi, identified 
equation (1.1) as presumably integrable in relation with the discrete Pm, which was 
called d-Pin (but should have been called (/-Pm for reasons that will become obvious). 
Indeed, as was shown in [4], the first form of the discrete Pm ever obtained: 

_ cd(x — a) (x — b) , . 

X ~ = (x-c)(x-d) ' { } 

with c, d constant and a, 6 proportional to A n , is associated to a linear problem of 
(/-difference type. 

In section 2 we will show the precise relation between (1.1) and (1.2). The method 
used for the derivation of q-Pyi is the "singularity confinement" discrete integrability 
criterion. In section 3 we present the bilinearization of both (/-Pm and q-Pyi equations: 
the former is obtained through the use of four r-functions while for the latter 8 r- 
functions are necessary. Section 4 is devoted to the Schlesinger transformations of both 
equations. 

2. Derivation of the q-Pyi equation from singularity confinement. 

We start with the g-Pm equation written in this general form (where a, 6, c, d, g may a 
priori all depend on n): 

_ g(x-a)(x-b) 
X ~~ ( x -c)(x-d) [ > 

In order to obtain the precise integrable form of this equation we use the singularity 
confinement approach. This method is based on the observation that, for an integrable 
mapping, the singularities that appear at some iteration, due to the particular initial 



conditions chosen, i.e. the 'movable' singularities, disappear after some iterations: they 
are "confined" . Based on this idea we can ask when a divergence appears in the iteration 
of (2.1). This happens whenever x passes through one of the roots of the denominator. 
Let us first consider the case where x = c while x is free. We find thus that x diverges 
and xc = g. Iterating further we would have found that the x's obtained do not in 
general depend on the free quantity x. The only way to restore this dependence, i.e. to 
confine the singularity, at this stage, is to balance the singularity of a; by a singularity 
of the rhs of (2.1). Thus x must be equal to a root of the denominator. The assumption 
x = c is not acceptable because it corresponds to a periodic singularity and is thus in 
contradiction with the requirement that the singularity be movable. We are left with 
the choice x = d and we have cd = g. In a symmetric way, starting from x = d we find 
cd = g. Thus the first condition for confinement is 

c/d = c/d (2.2) 

and c/d is periodic of period two. 

This is not, however, the only source of divergence: x may be equal to one of 
the roots of the numerator. This would lead to a diverging x unless the divergence is 
balanced by an appropriate vanishing factor. Let us start with the case x = a. This 
leads to x = 0. In order to obtain a finite x we must have x = b. (As previously, the 
confining condition x = a is not acceptable because it corresponds to a non-movable 
singularity). We find that we need ab = gab/cd = ab where the last equation is obtained 
from the b. Thus the second confinement condition writes: 

a/b = a/b (2.3) 

So a/b is also periodic of period two. In addition, using the expression of g above we 
also have 

abed = abed (2.4) 

We remark at this stage that the two conditions (2.2-3) relate n's of the same parity. 
Thus, the parameters with even and odd indices are only related by (2.4). 

If one introduces the simplifying assumption that the ratios c/d and a/b are not 
just periodic but strictly constant (as was done in [4]), one finds simply q-Pm. A 
dependent variable transformation x(n) — > x(n)c(n)/C allows one to replace c(n) by 
the constant C (which we will denote simply by c from now on). Since c/d is a constant, 
d is now also a constant. From (2.4) we find that a and b are both proportional to A n 
for some A. However, this simplification is unnecessary. Proceeding in full generality, 
we first use the dependent variable transformation (for even and odd indices separately) 
and impose c = c and thus d — d. We now find that (2.4) becomes abed = abed and the 
same equation at step n + 1 is abed = abed i.e. 

a/a = %/b (2.5) 



Combining the latter with (2.3) implies that the ratio a/a = b/b must be a parity 
independent constant. Calling this ratio q we find that both a and b are proportional to 
q n . We have in fact am = q n a e , a2 n +i = q n a (thereby breaking formally the symmetry 
between even and odd) and the analogous expression for b. Moreover, (2.4) now writes: 

qa e b e c d = a b c e d e (2.6) 

If q = A 2 and a a = Aa e , b a = A6 e , we recover the g-Pm case. As a last step we 
separate the even and odd x's calling, for instance, the odd ones y with the redefinition 
x(2n) — > x(n), x(2n + 1) — > y(n), and similarly for the (a, 6, c, d) which at odd n's will 
now be called (p, r, s,£). Because of the redefinition of n, what was previously called x 
is now just x. We find two coupled equations for x and y: 

st(x — a)(x — b) , n „ . 

= ^-P)^-^) (2.76) 

where c,d,s,t are constants and a,b,p,r are proportional to q n . The constraint (2.4) 
now becomes: 

prcd = qabst (2-8) 

We recognize immediately in (2.7), with the appropriate renaming of the coeffi- 
cients, the system (1.1) that defines the q-Pyi equation. Thus, q-Pyi is indeed contained 
in an embryonic state in q-Pm. 

3. Bilinear forms for q-Pm and q-Pyi. 

Let us start with the bilinearization of q-Pm. We expect x to be given by a rational 
expression involving the r-functions. The latter being entire functions, the divergences 
of x will be associated to the vanishing of r-fuctions which appear at the denominator 
of x. In order to proceed, we need the precise singularity structure of the equation [5]. 
From the previous section, we know that a singularity appears when x passes through a 
root of the denominator. Then x diverges at the next step and the singularity becomes 
confined at the subsequent step through the requirement that x passes through the 
other root. Two singularity patterns exist, {c, oo,d} and {d, oo,c}. The existence of 
these two singularity patterns suggests the introduction of two r— functions F, G: 

FG FG 
* = c(l + -^)=d(l + ^). (3.1) 

The singularity patterns described above correspond to the cases where either F or G 
passes through zero. Starting from this ansatz the bilinear form of g-Pm was given in 
[5]: 

cFG - dFG + (c - d)FG = (3.2a) 



cd = = — — 

■(cFG-dGF) + (c-a)(d-b)FG + c(d-b)FG + d(c-a)GF = (3.26) 



c — d — 

The first equation comes just from equating the two expressions of x. However, this is 
not the simplest form. In fact, as we have seen, two 'potential' singularities exist when 
x passes through a root of the numerator. The patterns are {a, 0, b} and {b, 0, a} and 
we remark that no singularity actually develops. This suggests the introduction of two 
more, auxiliary r-functions J, K in the following manner: 

„ FG, 7/1 FG JK . . 

x = c(l + T =) = d(l + =-) = m (3.3a) 

1 \, A JK^ 1,„ JK. FG , OOI . 

These equations in turn imply the q-Pm equation (1.2) for x. By equating the ex- 
pressions of x and 1/x we have four equations for the four r-functions and thus the 
bilinearization is trivially obtained: 

c(FG + FG) = d(FG + FG) = JK (3.4a) 

- (JK + JK) = T (JK + JK) = FG (3.46) 
a b 

In order to obtain the bilinear form for g-Pyi we start with the same remark as 
for the nonlinear variable: separating odd and even x's we introduced a new variable 
y. Similarly, the r-functions F and G (and also J, K if we are talking about (3.4)) are 
split into even and odd. Let us keep the old names for the even ones, and for the odd 
introduce the new r-functions M, N, P, Q. We look now for x and y: 

x = c(1 + 7f=> = d(1 + To ) = to (3 ' 5a) 

and the analogous expressions for 1/x, 1/y. However, the expressions are slightly more 
complicated if we leave full generality to c, a 1 , s and t. By rescaling either x or j/ it is 
always possible to set 

ca 1 = st, qab = pr, (3.6) 

which we will assume from now on. Rewriting (3.2a) for both even and odd indices we 
find: 

cMN - dMN + (c - d)FG = 

sFG - tFG +(s- t)MN = (3.7a) 

while (3.2b) becomes: 
cd — — 

' (sFG - tGF) + (c-a)(d- b)FG + c(d - b)MN + d(c - a)NM = 



s-t 



St 



(cMN - dNM) + (s-p)(t- r)MN + s(t - r)FG + t(s - p)FG = (3.76) 



c — d 



Similarly (3.4) can be transposed to the case of the 8 r-functions: 



c(FG + MN) 



d(FG + MN) = JK 



s(MN + FG) 



t(MN + FG) = PQ 

\(JK + PQ)=FG 
o 



1 



(JK + PQ) 



a 



1 



(PQ + JK) 



(PQ + JK) 



(3.8) 



MN 



P 



where (3.6) has been assumed. This completes the bilinearization of g-Pyi- 

4. MlURA AND SCHLESINGER TRANSFORMS FOR q-Pm AND g-Pyi- 

In this section we shall derive the Schlesinger transformations for q-Pm and g-Pvi- In 
the continuous case the Schlesinger are particular auto-Backlund transformations. As 
such they relate solutions of the same equation. However, the Schlesinger transforma- 
tions relate solutions corresponding to the same monodromy data except for integer 
differences in the monodromy exponents [2]. In the discrete case the relation of the 
Schlesinger transformations to monodromy exponents is not always clear. However, we 
can use an analogy with the continuous case. If one uses the proper parametrisation of 
the equation the Schlesinger transformations can be shown to be associated to simple 
changes of the parameters. The discrete case can be analysed in the same spirit. By 
using the proper parametrisation, one can identify, among the auto-Backlund transfor- 
mations, those which correspond to simple changes of the parameters and which can 
thus be dubbed Schlesinger's. As in the previous section, we shall start with g-Pm and 
then generalize the results to g-Pvi that is viewed as an asymmetric g-Pm. 

In order to find the Miura transform of g-Pm, we use the standard 'trick' intro- 
duced in [6]. Starting from the variable x expressed in terms of two r-functions F and 
G, we introduce a new variable u which can be expressed in terms of only one r-function, 
say G. We remark that (x — c) is proportional to F/F while (x — d) contains F/F. 
Upshifting the last object and multiplying by the first allows to get rid of F entirely. 
We thus find a first relation between x and u: 



where u is given by u = cdGG/GG. This is the first half of the Miura. Eliminating F 
between the two bilinear equations (3.2) leads to a (hexalinear) equation for G which can 
also be expressed in u. The latter would be the 'modified' equation of g-Pm ('modified' 
in the sense of the relation that exists between the Korteweg-de Vries equation and the 
modified Korteweg-de Vries equation). However, it is simpler to obtain the second half 
of the Miura and proceed from there. As we have shown in [6], the complement of the 



u = (x — c)(x — d) 



(4.1) 



Miura involves a rational expression which must be homographic in both u and u. We 
find readily: 

uu/cd — u — u + cd — ab , A _ 

% = ~ zz (4.2) 

—u/d — u/c + c + d — a — b 

Eliminating x and x between (4.1), (4.2) and its upshift leads to an equation for 
•u, u, u. This equation is, after a change of variables, a discrete form of q-Py equation 
(although not all the parameters of a q-Py are present). Introducing U = u — cd, we 
find: 

(UU - \ 2 abcd)(UU - oW) = «" g + UML + ME ± "Ml ± M (4.3) 
In order to define a different Miura we can introduce the quantity: 

w={\/x-l/a){l/x-\/b) (4.4) 

(we recall that b = b(n + 1) = Xb(n)) which depends only on the r-function K but not 
J. The second half of the Miura, analogous to (4.2), is: 

_ -aw/X - bwX + l/a+l/b-l/c-l/d 
abww + 1/ab — 1/cd — wX — w_/ X 

Again eliminating x leads to an equation for w. Introducing W = w — 1/ab, we find for 
this equation: 

(wW l ^ ( ww 1 g + v am )(^ + y ac )( w + wg + im 

{ X 2 abcd R — abed' XabW + 1 

(4.6) 

The quantity W = U/Xabcd satisfies an equation, obtained from (4.2), namely: 

tffiy L-)(WW 1 )- g + g + g + 1/A5C) g + 

X 2 abcd — a6c<i Aa6W + 1 

(4.7) 

Equation (4.7) has the same form as (4.6) provided one introduces the parameters 

a = a\f\~, b = b/VX, c = cVX, d = d/VX (4.8) 



We define 



and 



w = W + 1/Xdb = u/Xabcd (4.9) 



—aw/X — bwX + 1/a + lib — 1/c — 1/d ,.«„s 

x = — z^-rz — — - — ~ — — — (4.10) 

abww + 1/ab — 1/cd — wX — wj X 



Given this definition of x and since W satisfies (4.7) it follows that 



w = (i/x-i/a)(i/#-i/6) 



(4.11) 



and therefore x satisfies q-Pm with parameters a, 6, c, d. The transformation from x 
to x, through (4.1), (4.8), (4.9) and (4.10) defines an auto-Backlund transformation for 
g-Pin. In this case this is indeed a Schlesinger transformation which we denote by S%. 
(The convention used here is to give explicitly the parameters associated to x, rather 
than x, in (4.1) and (4.4).) The inverse transformation (5"") -1 can be obtained by 
defining w through (4.4), u = wXabcd and finally x through the analog of (4.2). 

In a similar way we can introduce the transformations S°[ = (S^) -1 and 
= (S 1 ") -1 . They correspond to multiplying the two parameters which appear ex- 
plicitly by \f\ while dividing the two others by the same quantity. These are the 
most elementary Schlesinger transformations. Using them we can construct further 
Schlesinger transformations that act separately on {a, 6} or {c,d}. For instance the 
product S®S2 corresponds to a — > aA, b — > 6/A, c — > c, d — > d. 

We turn now to q-Pyi and its Schlesinger transforms. Given its relation to g-Pm 
it is straightforward to find a Miura transformation analogous to (4.1). Note however 
that we need a two-component Miura: 

u = (x - c)(y - t), v = (y - s)(x - d) (4.12) 

and its second half: 

uv — sdu — ctv + st(cd — ab) , , _ „ . 

x = — = ^ -f (4.13a) 

— su — tv + st(c + d — a — b) 

uv — sdu — ctv + cdist — pr) 

V = — ; 77 ^ (4.136) 

—du — cv + cd(s + t — p — r) 

Eliminating x and y between (4.12) and (4.13) leads to an equation for u,v,v while 
eliminating x and y from (4.12), (4.13b) and the upshift of (4.13a) leads to an equation 
for u,v,u. We introduce, just as in the case of o-Pin, the new variables U = u — ct, 
V = v — sd and we obtain: 

/TTTr , , sd(U + cr)(U + bt)(U + cp)(U + at) 

(UV - abst) {UV - prcd) = — — — — (4.14a) 

U + ct 

^ rTT a/ta77 2i . ct(V + qsb)(V + rd)(V + qas)(V + pd) 

(VU - prcd) (VU - q 2 abst) = — ^— ^ — - — — (4.146) 

In parallel to (4.4) we introduce the Miura: 

w = (1/x - l/a)(l/y - 1/r), z = (1/y - l/p)(l/x - 1/6) (4.15) 

where now 6 = qb, the second half being: 

-qrw — pz + q(l/a + 1/6 — l/c— 1/d) 



x = 



prwz — qrw/b — pz/a + q(l/ab — 1/cd) 



= -aw - qbz + (1/p + 1/r - 1/s - 1/t) 
qabwz — aw/p — qbz/r + (1/pr — 1/st) 



The quantities W = w — 1/ar and Z = z — 1/ qpb satisfy the system: 



(WZ - M (WZ - l/ qabst ) = iW + V at) { W + y4r) { W^las ){ W + V a-) 
v — 1 /v ' ; pb(W + l/ar) 

(4.17a) 

, . , , . A (Z + l/pd)(Z + l/qU)(Z+l/pc)(Z + l/qbs) 
(ZW - l/ qabs t)(ZW - l/ qP rc d ) = 1 Jfz + l,J) — 

(4.176) 

In analogy to the g-Pm case we introduce W = U/Xabst and Z = V/Xprcd, (recall that 
q = A 2 ), and obtain: 

(WZ - q/prcdXWZ - 1/qabst) = *(W + X/dr)(W + 1/ Xbs)(W + X/pd)(W + 1/ Xas) 

/y 7 Aa6t(W + c/Aa6s) 

(4.18a) 



,~~ , w ~— , s t(Z + l/Xat)(Z + l/Xpc)(Z + l/Xbt)(Z + 1/Xrc) 

(ZW - l/qabst)(ZW - 1/qprcd) = -i ^ R UP± 1 R '- 

XprdyZ + s/ Xrpc) 

(4.186) 

This system has the same form as (4.17) provided one introduces the parameters: 



aps r 1 brc _ /pes 7 / adt . . 

° = »'— • 4= aV-- c = V"5"''' = VT ( ' 



/ape _ /6rs _ /acs ~ 1 /pdt 

P = A V~ r = V~ s = A VT' t= AVV 

We remind the reader that in section 3 we fixed the relative scaling of x and y so that 
cd = st and ga6 = pr (equation 3.6). It is easy to check that these relations are satisfied 
by the tilded parameters above. We define 

w = W + 1/ar = u/Xabst, z = Z + 1/qpb = v/ Xprcd (4.20) 

and x, y by the equivalent of (4.16) where w, z and all the parameters are tilded. Again, 
since W, Z satisfy (4.18) we have 

w = (l/x-l/a)(l/y-l/f), z = (1/y- l/p)(l/#- 1/6) (4.22) 

and x, y satisfy g-Pyi with 'tilded' parameters. The transformation from x, y to x, y is an 
auto-Backlund transformation B%£ for g-Pyi that is not a Schlesinger. (The convention 
again is that a are associated to x in w, p to y in z, c is associated to x in u and s to y in 
v). Starting from that auto-Backlund it is straightforward to obtain a Schlesinger 
transformation. In fact (B^) 2 is such a transformation where 

a — > aA, 6^6/A, c — > cA, d^d/A, p — > pA, r^r/A, s — > sA, £ — >£/A (4.23) 



However this is not the simplest Schlesinger transformation one can construct starting 
from B. In fact, through the product of auto-Backlund transformations B%[B%£ we 
obtain the parameters: 

/ abd } labc „ / bed ~ / acd 

p= aVV' r = A V~ s= aVV' f = 17 

Next we perform the transformation, X = y 1 abed/x, Y = ^prst/y and obtain q-Pyi 
for the variables A, Y. The resulting transformation is a Schlesinger corresponding to 
parameters: 

a — > a, 6^6, c — > c, d — > d, p — > pA, r^r/A, s — > sA, t— > t/A (4-25) 

Further Schlesinger transformations can be obtained through combinations of the ap- 
propriate B's. 

5. Conclusion 

In this paper we have undertaken a study of the discrete Pyi (q-Pyi) equation that 
was derived in [1]. We have shown that this equation can be viewed as a kind of 
'asymmetric' g-Pm. Indeed, its form is the most general one compatible with the 
singularity confienement requirement when applied to a g-Pm-type ansatz. Using this 
analogy we were able to obtain, in close parallel, results for the g-Pm and q-Py\. in 
particular we gave a simple bilinearization for both equations, involving four r-functions 
in the g-Pm and eight r-functions in the q-Pyi case. By deriving the appropriate Miura 
transformations for both equations we obtained their Schlesinger transformations. 

Several open questions remain at this stage. For example, the structure of the 
space of particular solutions of q-Pyi has only been touched upon in [1]. Rational 
solutions have not been considered at all. The coalescence chain of q-Pyi which is 
expected to contain new discrete Painleve equations has not been worked out yet. We 
hope to return to these questions in future publications. 
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